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Abstract 

We consider the bounds imposed by naturalness on the masses of superpart- 
ners for arbitrary points in nonminimal supersymmetric extensions of the 
standard model and for arbitrary messenger scales. We discuss appropriate 
measures of naturalness and the status of nonminimal supersymmetry in the 
light of recent experimental results. 



1 Introduction 



The problem of electroweak scale naturalness provides perhaps the most im- 
portant motivation for the consideration of supersymmetric extensions of 
the standard model (SSMs). Softly broken supersymmetry stabilizes the 
electroweak scale against quadratically divergent radiative corrections, but 
not against logarithmic divergences and finite corrections on the order of the 
superpartner masses. It is thus necessary that supersymmetry (SUSY) be 
effectively restored at a scale not much higher than the electroweak scale in 
order to avoid fine-tuning. 

Previous studies [pi B] have derived bounds on superpartner masses by 
requiring that there be no fine-tuning among the GUT-scale parameters of 
the supergravity-inspired minimal model with universal scalar and gaugino 
masses. The authors of relaxed this universality assumption slightly, al- 
lowing for nonuniversal scalar masses, but maintaining the assumption of 
gaugino mass universality and the assumption that soft SUSY-breaking pa- 
rameters are set at the GUT scale. In this paper, we consider the degree of 
fine tuning present at arbitrarily nonminimal points in the SSM parameter 
space and make no assumption regarding the so-called messenger scale at 
which soft SUSY-breaking parameters are set. 



1.1 Fine-tuning denned 

In the context of the SSM, the electroweak scale is determined at tree level 
by three mass parameters (m 2 ,, m 2 , m 2 ,) according to the equations 

m\ ml - ml tan 2 (3 

2 tan 2 /3-l 1 ' 

2m 2 7 , s 

sin 2/3 = z 2 (1.2) 

m x + m z y 

Expressions for these mass parameters in terms of SSM parameters will be 
given in section |]. Since the mass of the pseudoscalar Higgs boson m\ = 
m 2 x + m 2 , we can exchange the three mass parameters for m^, m^, and tan /3. 
Each point in (m z , mA, tan/3) space represents a possible Higgs sector of the 
SSM, and the actual Z boson mass defines a surface in this space. 

Qualitatively, the electroweak scale is fine-tuned if a fractionally small 
change in the value of a parameter results in a fractionally large change in 
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rriz- The degree of fine-tuning in a parameter a can thus be quantified by a 
sensitivity parameter A 

A = ^4f (1.3) 
m z oa 

which was first introduced by Barbieri and Giudice |1]. Bounding the allowed 
degree of fine-tuning, then, corresponds to requiring that the sensitivity A of 
mz to any parameter under consideration be less than some specified value; 
traditionally, one has taken A < 10, corresponding to 10% fine-tuning, as 
an acceptable level. Inserting equations ( |TTl~l - |1.2|) into (|1.3|) , we obtain the 
expression 
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We will impose the requirement A < 10 to place limits on the parameters of 
the SSM as functions of the messenger scale, tan/?, and tua- 

As pointed out by Dimopoulos and Giudice @, the fact that the coeffi- 
cients of 5m 2 , and 8m 2 vanish in the tan (3 — > oo limit implies that m y must 
exhibit a dependence on a parameter in order for the requirement of a partic- 
ular degree of naturalness to impose a globally valid limit on the parameter. 
This situation arises because the decoupling limit 

m x ~mz tan ft m y ~ mz m z ~ m z tan 1 / 2 f3 (1.5) 

is a perfectly natural corner of parameter space with large tan f3. The form 
of the dependence of m x and m z on a parameter are nonetheless relevant in 
determining how a naturalness bound varies with tan (3 and m^. 

We should emphasize that the degree of fine-tuning present in a particu- 
lar realization of the SSM depends upon the choice of a parameter set which 
we regard as specifying the theory. Previous studies have taken as a pa- 
rameter set the DR soft-SUSY breaking parameters at the GUT scale. The 
messengers which communicate SUSY-breaking to the SM particles and their 
superpartners could in principle, however, exist at any scale. Whatever the 
scale at which the SUSY-breaking parameters are set, their effect on the elec- 
troweak scale is determined by integrating the renormalization group (RG) 
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equations which describe their flow from the messenger scale rh down to the 
electro weak scale m z ^\ 

In section |2], we integrate the relevant one-loop RG equations in order 
to find the dependence of the electroweak scale on the the DR parameters 
at the messenger scale, independent of any assumptions of universality. In 
section [3], we apply equation (|1.4j ) to obtain limits these on parameters and, 
in turn, on sparticle masses as functions of the messenger scale. Finally, in 
section [|, we consider mass limits on sparticles which first contribute to rriz 
at the two-loop level. 

1.2 Naturalness vs. sensitivity 

The equation of sensitivity with unnaturalness has been criticized by An- 
derson and Castano Pj, who rightly point out that it is inappropriate in 
some cases. For example, a small scale m = Ae _47r / 9 resulting from dynam- 
ical symmetry breaking exhibits a strong sensitivity to g even for perfectly 
natural values of g ~ 1. This motivates them to introduce a more refined 
measure of naturalness which, instead of simply reflecting the sensitivity of 
mz to an underlying parameter, compares that sensitivity to the average 
sensitivity over the entire allowed parameter space. We do not employ their 
refined definition for this study. We do, however, present in this section some 
considerations which are, to our knowledge, original and which we believe 
elucidate the relationship between naturalness and sensitivity and clarify the 
circumstances under which their equation is appropriate. 

Suppose an underlying parameter set a determines the electroweak scale 
according to a function m(a) and we imagine that the underlying parame- 
ters are distributed according to a probability distribution p(a). Then the 
probability distribution of electroweak scales is given by 

p(m') — J dap(a) 5(m' — m(a)) (1.6) 

and the probability of the electroweak scale being as low or lower than the 

1 In fact, this procedure only yields the contributions to mz enhanced by factors of 
\a{m/mz). Unenhanced terms, while present, are not only presumably smaller, but also 
scheme-dependent. We therefore ignore the unenhanced terms, but note that the the 
limits derived here may therefore be considered valid only for messenger scales satisfying 
\n(rh/mz) ^S> 1. 
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measured value m is 

P(m) = / dm' p(m') (1-7) 
Jo 

We typically imagine the probability density of parameters p(a) to be some 
flat function over an allowed range. Given some p(a), the probability P(m) 
would seem an appropriate measure of the naturalness of a low electroweak 
scale. 

In a multidimensional parameter space like that of the SSM, the exami- 



nation of the P{m) implied by equations ([176 -1.7) for various p(a) seems a 
daunting task. We can nevertheless glean some insight by considering the 
simplified case in which m depends depends only on a single parameter a, 
distributed uniformly between between and a max . Then 

/ i \ 1 / 5m\ 1 
P(m') = — (1.8) 

Amax V Oa / 0=m -l( m /) 

and, assuming 5m/ 5a to be relatively constant over the allowed values of a, 
we obtain 

m 5a 



P m) 



a max 5m 



;i.9) 



which implies A ~ 1/P, i. e. sensitivity is the inverse of naturalness. In 
the example of Anderson and Castano, it is the strong dependence of 5m/ 5g 
on g which violates the assumption made for our toy example and belies 
the simple relationship between A and P. Moreover, the integrations de- 
scribed by equations ( |1.6| - |1.7| ) may be easily carried out for the Anderson 



and Castano's example to find that, for g distributed uniformly on the in- 
terval < g < 0(1), a hierarchy of scales m/A ~ 10~ 10 is not terribly 
improbable, while hierarchy m/A ~ 10~ 100 is quite improbable. 

We believe that equations ( |L6| - |i.7| ) provide a quantitative definition of 
naturalness similar in spirit to that of Anderson and Castano, but concep- 
tually simpler, since it allows naturalness to be defined a priori without 
reference to sensitivity. Our experience with a toy example has suggested 
that the qualitative equation of sensitivity to unnaturalness is justified when 
the value of 5m/ 5a is not strongly dependent on the location in parameter 
space a. Since such an approximation is indeed valid for the parameters 
considered here, we believe our equation of sensitivity with naturalness to be 
qualitatively valid for the case under consideration. 
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1.3 The calculation sketched 



Before embarking on this calculation, we outline a didactic method which 
yields back-of-the envelope estimates of the limits set by naturalness and 
the form of their dependence on A and the messenger scale m. If a mass 
parameter m contributes to a higgs mass at tree level, then we expect from 
naturalness that. 

0(l)m 2 <m 2 z A (1.10) 
m < O (300 GeV) J (1.11) 

If, on the other hand, a particle with mass m contributes at one-loop via 
a coupling g ~ 1, we expect 

oMjfvi^i-) ™ 2 ;$™I A (i-i2) 



{An) 2 \m 



/33 A \ 1/2 
m < 0(650 GeV) (j- J (1.13) 



where t = ln(fh/mz), so that t ~ 33 for a GUT messenger scale. As this 
calculation indicates, unlike the bounds on parameters entering at tree-level, 
bounds on the masses of particles contributing at one loop depend on the 
messenger scale. 

Finally, if the particle contributes at two loops, we expect 

° (1) (l^ ln (l) m2 - m|A (L14) 

/ 33 A \ 1/2 
m<(8000GeV)(^ T -J (1.15) 

which is, apparently, a significantly weaker bound. 

From these considerations, we expect bounds on mass parameters which 
enter into higgs masses at tree level to be messenger scale independent while 
those which enter as radiative corrections to depend on the messenger scale 
approximately as ln^ 1 ^ 2 [mjmz). All mass bounds are expected to scale as 
A 1 / 2 . Our work in the following sections will bear out these qualitative 
expectations. 
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2 One-loop RG equations and their solutions 



The tree-level relations between the parameters (m^,m^,m^) and the soft 
super symmetry breaking Higgs masses and higgsino mass parameter /i are: 

ml = m 2 hd +^ (2.1) 
m 2 y = ml u +f (2.2) 
m 2 z = m 2 ud (2.3) 

We are thus required to integrate the RG equations || for four mass param- 
eters — rrih d , rrih u , fi and m U( i — in order to determine their dependence on 
messenger scale parameters. The particles coupled to the higgses by gauge 
interactions or the top Yukawa coupling at one loop are: gauginos, higgsinos, 
left- and right-handed stops and left-handed sbottomsQ. Since these couplings 
are large, we expect the masses of these particles to strongly renormalize the 
higgs masses; examination of the relevant one-loop RG equations bears out 
this expectation. Because of their strong effect on the higgs masses, natural- 
ness places strong constraints on the masses of these particles; for this reason, 
the authors of @ referred to these particles as "brothers of the Higgs" . 

We begin with the equations for rrif ld and m^. The RG equation govern- 
ing the evolution of the down-type Higgs' mass squared reads 

odm 2 , 

( 4 *) ~dT = ~ fh * {2A) 

while the up-type Higgs mass squared evolves according to the coupled equa- 
tions 



(2.5) 





m lu\ 




f fh u 


i , 


l-l 


u 


™ 2 U J 







2 For tan/? > mt/rrib, the bottom quark Yukawa coupling and potentially even the 
tau Yukawa coupling are large enough to require consideration. Since this represents an 
extreme corner of parameter space and would much complicate the RG equations, we do 
not consider this possibility here. We also ignore the effects of soft-SUSY breaking trilinear 
couplings, which may be present in principle, but are small in most models. 
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Here the inhomogeneous terms are 

6 



fh d 
fh u 

fu 



6g 2 2 m 2 2 + -glm\ + -g\M 2 



r 2 2 i 2 2 
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The gauge couplings and gaugino masses evolve according to 



9i(t) = 9aZi{t) m{t) = {mi) z i{t) 



where 



Zi{t) 



(4tt) 2 ' 



-Bl,2,3 



33 



(2.6) 
(2.7) 
(2.8) 
(2.9) 

(2.10) 
(2.11) 



The mass M, defined by 

M 2 = 2 tr [Fm 2 ] = m 2 u - m 2 hi + tr [m 2 - m 2 - 2m 2 u + m 2 d + m 2 ] (2.12) 



may be taken to be a RG invariant, since the terms proportional to gaugino 
masses vanish in its evolution equation, while the presumably dominant first 
and second generation scalar masses are only weakly renormalized by Yukawa 
couplings. 

Before continuing, let us clarify our notation. Parameters without any 
additional annotation (e. g. g or m) refer to parameters evaluated at any 
scale t = ln(m/mo). Parameters annotated with a twiddle (e. g. m) refer to 
the value of the parameter at the messenger scale t = ln(m/mo), and thus 
constitute what we will regard as the SSM parameters a in the context of 
equations ( |1.3| ) and ( |1.4j ). Finally, parameters with a nought (e. g. go) refer 
to GUT-scale parameters evaluated at t = 0. 

~4|) can be readily integrated 



The equation for the evolution of m 2 
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(2.13) 
(2.14) 
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The solution of the coupled equations ( |2.5| ) governing the evolution of m hu is 
more complicated but can also be written in closed form in a useful approx- 
imation outlined in an appendix. 

The RG equation governing the evolution of the /i-parameter, 



2 dfi 



may be straightforwardly integrated to yield 



(2.15) 



(2.16) 



Here g(i) is a function of gauge and Yukawa couplings whose specific form is 
irrelevant for our purposes. 



Finally, the RG equation for the evolution of rn^ d , 



(4tt) 



dm 2 
~dt 



ud 



[2.17) 



may be likewise straightforwardly integrated 



m 2 s + I 3 : 



m 2 + 



\Z X -Z X „ 



11 z x 



-mi fi 



[2.1%) 



3 Limits on brothers' masses 

Having solved the RG equations for the relevant parameters, we can proceed 
to apply equation (|1.4j) to derive bounds from naturalness on the parameters 
which enter into their solutions. 



3.1 Gluino masses 



The algebra involved in limiting gluino mass parameter 7713 is easiest, because 
it enters only into the expression ( |A.2| ) for rrih u . Applying the naturalness 
criterion (|1.4j) yields 



32 ml 
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9 m 2 z Z\ 
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Figure 1: A = 10 limit on gluino mass as a function of the messenger scale 
in the large tan/3 limit. The scaling of this limit with tan/3 and itla is given 
by equation fl3.1|) . 



where c = 1 + m 2 z /m\ > 1, and Zj, A and are defined by equations 
(|2.11| ), ( |A.6|) , and ( |A.7|) , respectively. We have used equation (|2.10|) to con- 



vert a limit on m 3 to one on m 3 itself. Note that the limit becomes signifi- 
cantly more stringent for smaller values of tan (3. Figure |l] shows the variation 
of the limit with the messenger scale in. 

Because current experiments |§ already require > 200 GeV in the 
context of the minimal model with universal GUT-scale scalar and gaugino 
masses, our work allows us to deduce that the gaugino mass in the MSSM 
is fine-tuned to at least 11%; this required degree of fine-tuning increases by 
more than a factor of two as tan/3 is reduced to tan/3 ~ 2. 

3.2 Chargino and neutralino masses 

The analysis of the limits on the bino, wino, and higgsino mass parameters 
mi, m 2 and (jl is considerably complicated by the fact that each of these 
parameters enters into more than one RG equation, and in such a way that 
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the application of equation (O) results in inequalities involving all three, so 
that we must solve three inequalities simultaneously for each value of tan /3 
and rriA in order to derive bounds on the parameters. 

While this procedure can in principle be carried out numerically, the 
variation of the resulting bounds with tan/3, tua, and the messenger scale in 
cannot be displayed in a compact form. We therefore content ourselves with 
the presentation of bounds in the large tan/? limit, in which the relevant 
inequalities decouple, and with the observation that, as was the case for 
gluino masses, mass bounds generically become more stringent as tan /3 is 
lowered. 

In the large tan/3 limit, as noted previously, only the dependence of w? y 
on the parameters matters, and the following inequalities are readily derived: 
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(3.2) 
(3.3) 
(3.4) 



These limits are displayed in figure |2| as a function of messenger scale. Most 
interesting is the limit on u, first because it is the most stringent and sec- 
ond because, as we expected for parameters which enter at tree level, it is 
independent of the messenger scale. 

Unlike the gluino mass m^, the parameters mi, mi, and a are not them- 
selves particle masses, but enter into mass matrices which must be diagonal- 
ized in order to determine particle masses. As pointed out by Barbieri and 
Giudice |T|, however, from the structure of the relevant mass matrices it fol- 
lows that the masses of the lightest chargino and neutralino are constrained 
by 



m x ± < m w + mm (m 2 , a ) 
mt 



^ ■ ( 2 i 2-2/1 2 i 2 2 q 2,1 2 \ 

o < mm (mi + m z sin tV, m 2 + m z cos uw,u + 2 m z) 



(3.5) 
(3.6) 



Using these inequalities, our upper bounds on these parameters may be trans- 
lated into upper bounds on the masses of the lightest chargino and neutralino. 
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Figure 2: A = 10 limits on wino and higgsino mass parameters as a function 
of the messenger scale in the large tan/3 limit. The corresponding limit on 
the bino mass parameter is greater than 1 TeV. 
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In all cases, it is the bound on // which dominates, implying 

m x ±<165GeV (3.7) 
m x o<160GeV (3.8) 

Because they arise from the bound on /j, these bounds are independent of 
the messenger scale. Although they do not scale exactly as A 1 / 2 , they do so 
to a good approximation for A > 10. 



3.3 Scalar masses 

Now consider the limits on the scalar mass parameters rhh u , rh q , and fh u . 
From their contributions to expression ( A.2 ) for rrih u and the application of 
the naturalness criterion ( |1.4| ), we obtain 
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The physical masses of the left- and right-handed stops, as determined 
by equations (|A.3|) and ( |A.4| ), receive contributions from both scalar and 
gaugino masses. In the case where the contributions from scalars dominate, 
it is easy to translate our bounds on scalar mass parameters at the messenger 
scale into bounds on the stop masses. 
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Figure 3: A = 10 limits on left- and right-handed stop masses in the large 
tan/? region as a function of the messenger scale. The solid lines show the 
limits implied when scalar masses are large compared to gaugino masses; 
when gaugino masses are comparable, the dashed lines result. The scaling 
of the solid lines with tan/3 and tua is given by equations ( p.ll| - |3~T^) . The 



bump in the mass limit on right-handed stops near 10 13 GeV is due to a 
positive contribution from the bino mass parameter, which is allowed to be 
large due to a zero in the coefficient with which it contributes to m hu near 
this scale. 

These bounds are shown as solid lines in figure ^. It is noteworthy that both 
of these bounds cannot be saturated simultaneously, since m 2 contributes 
negatively to m? q and vice versa. The scaling of these limits with tan /3 makes 
them more stringent by more than a factor of two for low values of tan (5 ~ 2. 

When gaugino masses are comparable to stop masses, their possible con- 
tributions must also be considered. Inserting the gaugino mass limits ( |3 . 1| ) 
and (|3.2| - |3.3|) into equations (|A.3|) and (|A.4j ) gives the stop mass limits indi- 



cated by the dashed lines in figure ||. 

The scalar mass sum M 2 enters into the expressions ( p. 14 , |A.2| ) for both 
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Figure 4: A = 10 limit on generic sparticle masses in M ^ models in the 
large tan 3 regime as a function of the messenger scale. The scaling of this 
limit with tan 3 is determined by equation ( [3. 13] ). The limit for right-handed 
up-type squarks is more stringent by a factor \[2. 

rrih d and vnh u . Treating it as a parameter gives a naturalness limit 

I M 2 , ( Z x \ tan 2 3-1 ,„ „. 
— —5- In - 1 <A ^ — 7 3.13 

II m| I Zi / tan 2 /3+l v ; 

on M and by implication on the masses of the scalars which enter into it. 
The variation of this limit with the messenger scale m is shown in figure f| 
in the large tan/3 limit; again, the limit becomes more stringent as tan/3 
becomes smaller. 

Of course, M itself should not be regarded as a parameter but as a func- 
tion of the scalar sparticle masses. As noted by Dimopoulos and Giudice , 
if some symmetry, e. g. scalar mass degeneracy or SU (5), ensures M = then 
equation ( |3.13| ) gives no limit on any scalar masses. If, on the other hand, 



scalar masses are not so constrained, this is an upper limit on the masses of 
the scalar partners of fermions of all three generations. 
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4 Cousins of the Higgs 

If some symmetry principle does protect M = 0, the masses of the first and 
second generation sparticles enter the one-loop RG equations suppressed by 
tiny Yukawa couplings and their masses are thus essentially unlimited by one 
one-loop radiative corrections. On the other hand, at two loops the inhomo- 
geneous terms in the RG equations described above get contributions 



3 \9tSl + 3g$Si - \g\T 2 (4-1) 

(47r)4 ^r D l atsl + 3 ^ 2 + l alT2 (4 - 2) 

(47r)4 ^r D h gts! + 3glsi + t 9 * 81 + \ 9lT2 (4 - 3) 

(4vr) 4 ^ D ±gt£% + fgtS! - jj«#Z* (4.4) 

where 

S? = 2§tr [U(1)W] 

= I (3m 2 hu + 3m 2 hd + tr [m 2 + 8m 2 u + 2m 2 d + 3m 2 + 6m 2 e ] ) (4.5) 

S 2 2 = 2 tr [SU(2) 2 m 2 ] = m\ u + m 2 hd + tr [3m 2 q + mf ] (4.6) 

S 2 = 2 tr [SU{3) 2 m 2 ] = 2m 2 q + m 2 u + m 2 d (4.7) 

and 

T 2 = g\Tl + g \T 2 2 + g\Tl (4.8) 

T i = To « ~ O + tr lh m l ~ iK + S™2 ~ To m i + l m l] ( 4 -9) 

T l = I «, - m 2 hd + tr [m\ - m 2 ]) (4.10) 

T 2 = ftr [m 2 - 2m 2 u + m 2 d ] (4.11) 

Including these mass sums gives a contribution to the up-type higgs mass 

m 2 hu D |° (C.S 2 + C 2 S 2 2 + C 3 S 3 2 + D,T 2 + D 2 T 2 + D 3 T 2 ) (4.12) 
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where 
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(4.13) 
(4.14) 
(4.15) 

(4.16) 

(4.17) 
(4.18) 

(4.19) 



Equation ( 4.12| ) implies naturalness limits for all sparticles of all generations 
given by 
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\C S \ < A 



where the coefficients C s may easily be read off as 

C q = id + 3C 2 + 2C 3 + i£>! + %D 2 + |D 3 
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(4.20) 



(4.21) 
(4.22) 
(4.23) 
(4.24) 
(4.25) 



The resulting limits are shown in figure [5[ This shows that it is possible 
to raise the masses of the q- and d-type first and second generation squarks 
to between 4 and 8 TeV, and those of the other squarks and sleptons even 
higher, without violating naturalness. 
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Figure 5: A = 10 limits on sparticle masses from two-loop naturalness in the 
large tan (3 limit as a function of the messenger scale. 
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5 Conclusions 



We have derived bounds on the masses of supersymmetric partners from 
the requirement of naturalness in arbitrarily nonminimal incarnations of the 
SSM in which SUSY-breaking is communicated by messengers at an arbitrary 
scale. 

Examination of the variation of naturalness bounds with messenger scale 
reveals that these bounds may be increased from their GUT messenger val- 
ues by between 15% and 100% by lowering the messenger scale for all super- 
partners except the lightest charginos and neutralinos, for which we obtain 
messenger scale independent mass bounds of ~ 160 TeV for 10% fine-tuning. 

We have derived the variation of several of these several mass bounds 
with tan (3 and found that naturalness more significantly constrains sparticle 
masses in models with low tan (3 than in models with high tan j3. The mass 
constraints in the two regimes can differ by more than a factor of two. 

The most problematic constraint for the traditional MSSM with GUT 
scale mediation of SUSY breaking is that placed on the gluino mass m 3 < 
260 GeV. Present experiments thus require that the MSSM be at least 11% 
fine-tuned and more than 5% fine-tuned for tan/5 < 2. This problem can be 
evaded by models which significantly lower the messenger scale or make the 
gluino the lightest supersymmetric particle. 

As accelerator searches progress, this work will place increasingly strin- 
gent constraints on possible realizations of the SSM. 
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A Appendix 



Equations (|2.5| ) represent an inhomogeneous linear differential system which 
may be solved formally in closed form. In order to make computational sense 
of this formal solution, we must examine the behavior of the top Yukawa 
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coupling A(t), which evolves according to the RG equation 
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This equation lacks a closed-form analytic solution, but an nonetheless be 
integrated numerically. Numerical integration indicates that the top Yukawa 
coupling remains within ten percent of A ~ 1 over nearly^ the entire in- 
teresting range of t and tan/?. We will therefore perform the integrations 
considering A to be constant. This yields the results 
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3 The exception occurs for GUT scale values of t when tan/3 < 1.5. Even an exact 
treatment of the evolution of A would not be particularly helpful in this case, since its 
behavior in this region becomes highly sensitive to the exact value of nit- 
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where 
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